The friction factor of two-dimensional rough-boundary turbulent soap film flows 
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We use momentum transfer arguments to predict the friction factor / in two-dimensional turbulent 
soap-film flows with rough boundaries (an analogue of three-dimensional pipe flow) as a function 
of Reynolds number Re and roughness r, considering separately the inverse energy cascade and the 
forward enstrophy cascade. At intermediate Re, we predict a Blasius-like friction factor scaling of 
/ oc Re"^''^ in flows dominated by the enstrophy cascade, distinct from the energy cascade scaling 
of Re~^/*. For large Re, / ~ r in the enstrophy-dominated case. We use conformal map techniques 
to perform direct numerical simulations that are in satisfactory agreement with theory, and exhibit 
data collapse scaling of roughness-induced criticality, previously shown to arise in the 3D pipe data 
of Nikuradse. 
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Turbulent flows are marked by rich structure over a 
range of scales — they host fluctuations, vortices, tangles, 
and other coherent structures that continue to defy a 
detailed, analytical understanding [l|, When parame- 
terized in terms of the typical flow speed U , character- 
istic length scale L and kinematic viscosity of the fluid 
V, three-dimensional turbulence exhibits universal phe- 
nomena as the Reynolds number Re = UL/v — > oo. 
Most famously, in a theory referred to as K4l[^ the 
dependence of the fluctuation energy spectrum E{k) on 
wavenumber and mean energy transfer rate e occurs in 
a way that is independent of v: E{k) — £2/3^-5/3 
values of wavenumber in the so-called inertial range, in- 
termediate between the scales of forcing and the scales 
where molecular viscosity becomes significant. In this in- 
ertial range, turbulent eddies break up into smaller eddies 
through a mechanism which is to a first approximation 
Hamiltonian, and results in a cascade of energy to smaller 
length scales |5j. 

During the 1930's, Nikuradse undertook a systematic 
series of measurements of the pressure drop across a tur- 
bulent pipe flow as a function of ReQ and also as a func- 
tion of r/R, the scale of the roughness of the pipe walls 
r, normalized by the pipe radius i?[7]. The former mea- 
surements provided strong support for Prandtl's bound- 
ary layer concept Q, and have been replicated and sur- 
passed only recently[9|, while the latter measurements, 
despite recent efforts|lCllll|. remain to this day the most 
complete data set of its kind, spanning three orders of 
magnitude in Reynolds number and a decade in the di- 
mensionless roughness r. These data reveal that the fric- 
tional drag experienced by a turbulent fluid in a pipe with 
rough walls is a non-monotonic and complicated function 
of Reynolds number and roughness, which despite intense 
interest and practical importance (see, e.g. Ref. ^), has 
only begun to be understood 1^, ll3| through two related 
developments. 

First, Gioia and Chakrabortvfl^ estimated the 
momentum-transfer between the walls of the pipe and 



the flow, explicitly taking into account the presence of 
roughness. Their resultant formula for the dimensionless 
friction factor (defined precisely below) is expressed in 
terms of the turbulent kinetic energy spectrum E{k), and 
thus makes a direct connection between a macroscopic 
flow property and the velocity field correlations. Second, 
Goldenfeld [13| pointed out that the power law behavior 
of Nikuradse's friction factor data in the regimes Re — > oo 
and r/R — > was analogous to critical phenomena, where 
the inverse Reynolds number and roughness play similar 
roles to, for example, the coupling constant and exter- 
nal magnetic field in an Ising model. Consequently, the 
dependence of Nikuradse's data on Re and r can be col- 
lapsed onto a universal function with sufficient precision 
for intermittency corrections to be extracted IJ] . These 
results show that the friction factor reflects the nature of 
the turbulent state through its dependence on the energy 
spectrum, and that the turbulent state is itself a mani- 
festation of a non-equilibrium critical point at Re — oo 
and r/R 0. 

The purpose of this Letter is to test the claims of Refs. 
and [13] in a context where detailed calculations are 
in principle possible: the case of two-dimensional soap- 
film turbulence [isl. [l6| . Here, a soap film is supported 
between two vertical wires, and the draining flow pro- 
vides a versatile laboratory for exploring two-dimensional 
turbulence [Til . It is well- understood that the nature of 
turbulence in 2D is different from 3D: there is no vor- 
tex stretching, for example. Nevertheless, turbulent phe- 
nomena exist, and possess the novelty that there are two 
cascades: an energy inverse cascade that runs from small 
to large scales [171 ll8|, and a forward cascade[18] in the 
enstrophy ^l = |V x vp, where v is the fluid velocity 
field. This enstrophy cascade yields an energy spectrum 
E{k) ~ 0^/^k~'^, where (3 is the rate of transfer of en- 
strophy. 

Prior work, dating back to Prandtl and others (for a 
review see Ref. , Sj ) is not able to make a prediction about 
the friction factor in these cases, because it has no spe- 



2 



cific representation of the nature of the turbulent state, 
and in particular is disconnected from the energy spec- 
trum. On the other hand, the momentum-transfer theory 
of Gioia and Chakraborty[l3| can reflect the character 
of 2D turbulent states, as expressed by the energy spec- 
trum, through the dependence of the friction factor on Re 
and r. We show below that the momentum-transfer the- 
ory predicts a significant dependence of the friction factor 
on the nature of the turbulent cascade, one that we ob- 
serve in direct numerical simulations reported here, and 
which obeys the scaling predicted by roughness-induced 
criticality. Thus our direct numerical calculations agree 
well with the momentum-transfer and roughness-induced 
criticality picture, and strongly suggest that the standard 
picture of turbulent boundary layers is incomplete. 

Calculation of the friction factor scaling laws in 
2D:- In the momentum-transfer theory of Gioia and 
Chakraborty, the friction factor is shown to be propor- 
tional to the root-mean-square velocity fluctuation Us at 
a scale s determined by the larger of the roughness r 
or the Kolmogorov scale rjK- Since E{k)dk represents 
the turbulent kinetic energy in the wavenumber band be- 
tween k and k + dk, it follows that: 



E{k)dk 



l/s 



1/2 



(1) 



Anisotropy near the wall has only a small effect [19] on 
the low-order structure function used in our calculation. 
For simplicity, using the K41 form for Eik), we obtain 
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With the limiting forms for s at large 



and small Reynolds number, we obtain the predictions of 
the empirically-observed Blasius regime "2^, in which the 
friction scales as Re^^^^ for small but turbulent Reynolds 
numbers, and the Strickler regime 2]J at large enough 
Re, where the friction factor is independent of Re and 
only depends on the roughness through the relation / oc 
(r/i?)i/3. 

In two-dimensional turbulent systems, both the energy 
cascade or the enstrophy cascade may be observed, or 
they may occur individually [23| depending on the man- 
ner of energy injection and the scale at which it occurs. 
The two dimensional inverse cascade friction factor is 
the same as the case of three dimensional flows, with 
a Blasius scaling of / oc Re~^^^ and a Strickler scaling 
of / oc (r/Ry/^. The energy spectrum due to the en- 
strophy cascade leads to a new prediction for the friction 
factor: a scaling of / oc Re ' in the Blasius regime 
and / oc (r/R) in the Strickler regime. These are our 
central predictions, which we seek to verify by numerical 
simulation in the next section. In general, the friction 
factor corresponding to any conserved quantity (such as 
helicity) with units [L]"[T]'' is / oc Rc^^^^'^^/^^-'^) (Bla- 
sius regime) and / oc [r/R)^^'^ (Strickler regime), where 
</. = a/(l-6). 
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FIG. 1: Energy spectra for grid- and roughness-generated 
turbulence. Grid-generated turbulence exhibits the k~'^ en- 
strophy cascade, whereas roughness-generated turbulence ex- 
hibits the k~^^^ inverse cascade scaling. Inset: simulated wall 
velocity profile of grid-generated turbulence in a smooth pipe 
at Re = 60000. The profile is consistent with a power law 
with exponent 0.323 ± 0.005. We predict an exponent of 1/3 
for enstrophy cascade turbulence. 



Simulations of 2D turbulent rough-pipe flows:- To test 
the momentum transfer theory's prediction of the fric- 
tion factor in 2D, we have performed simulations for a 
range of Reynolds numbers and single- wavelength rough- 
ness, both with grid-generated turbulence and turbulence 
generated by wall roughness. The roughness of the wall 
breaks translational invariance and means that one can- 
not simply solve the Navier-Stokes equations using spec- 
tral methods. We have overcome this difficulty by using a 
judiciously-chosen conformal map technique, allowing us 
to use a spectral method to satisfy incompressibility. The 
SMART algorithm (2^ is used to calculate the advection 
of the velocity field. The friction factor is measured by 
computing the pressure drop necessary to maintain the 
average flow velocity over the periodic domain. 

To simulate a rough- walled pipe, we apply a conformal 
map of the form w = z + r eKp{ikz), where the aspect ra- 
tio is held constant (rfc = 3/4) and the wavenumber k 
may be varied to produce roughness of different scales. 
Note that r plays the role of roughness in Nikuradse's ex- 
periments, but our aspect ratio is 3/4 and not unity as in 
his experiments. This conformal map results in the addi- 
tion of two body force terms to the Navier-Stokes equa- 
tion in the transformed (rectangular) domain, in addition 
to an overall weighting factor deriving from the changed 
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FIG. 2: Scaling of the friction factor with respect to Re for 
inverse cascade and enstrophy cascade dominated flows in 2D. 
The roughness is r / R = 0.067, and the data have been aver- 
aged over a time of 5 pipe transits. 
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velocity in the transformed coordinates, and the vector 
A is defined as: 



(3) 



We use a simulation domain of 2048 x 512 to simulate 
a section of pipe of diameter 1 and length 4. After initial- 
izing the velocity field we allow the system to evolve for 
a sufficient number of pipe transits so that the system 
is fully turbulent (one pipe transit corresponds to four 
units of time as the mean flow velocity is set to 1 in the 
simulation units). The smaller the roughness, the more 
transits are needed. This results in roughness-generated 
turbulence, in which case the observed energy spectrum 
is dominated by the inverse cascade, as shown in Fig. ([1]) . 

In order to attain an enstrophy-dominated flow, we 
used a technique suggested by the observations reported 
by Rutgers |22|. We simulated grid-generated turbulence, 
by placing a series of cylinders at the mouth of the pipe; 
in each cylinder we set the velocity field to zero every 
timestep. After one pipe transit the velocity field is fully 
developed. We then remove the grid and allow the tur- 
bulence to decay for a transit before we begin to measure 
the friction factor and other fiow properties. We have ob- 
served energy spectra dominated by the enstrophy cas- 
cade in this system, as shown in Fig. ([T]). 

Our simulation results at small values of the dimen- 
sionless roughness (r/R = 0.067) are plotted in Fig. 
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FIG. 3: (Color online). The bottom inset shows the enstrophy 
cascade data collapse of the friction factor curves for nondi- 
mensional roughness 0.05 (o), 0.08 (□), 0.1 (A), and 0.2 (v) 
over a range of Reynolds numbers from 1000 to 80000. The 
top left insent shows the unsealed friction factor data. The 
top right inset shows the energy spectrum at r/R = 0.08 and 
Re = 80000. The straight lines correspond to k~^^^ and k~^. 



These results were obtained by averaging over 5 full 
pipe transits, yielding reproducible values for the fric- 
tion factor, with controlled error bars, as shown. For 
this flow we observe an approximate power-law scaling 
of the friction factor with Reynolds number, with an 
exponent —0.22 ± 0.03 together with an energy spec- 
trum dominated by the inverse-cascade. In the case of 
grid-generated decaying turbulence, corresponding to an 
enstrophy-cascade dominated spectrum, we observe an 
exponent of —0.42 ± 0.05. These results are within sat- 
isfactory agreement with the scalings of —1/4 and —1/2 
respectively, predicted for the 2D Blasius regime on the 
basis of a momentum transfer argument. 

We cannot reach sufficiently high Reynolds numbers 
to observe a pure Strickler regime, but we can verify 
the Strickler scaling exponent with data collapse. In 
three dimensions, or in a system dominated by the in- 
verse cascade, we expect data collapse when plotting 
/Re^/^ against {r/R)Re^^^^. For the enstrophy cas- 
cade, these variables should be /Re^^^ and (r/i?)Re^^^ 
respectively. We have observed previously that in the 
presence of roughness, the spectrum is dominated by the 
inverse cascade. However, we have found that by adding 
a small amount of random forcing to the velocity field, the 
enstrophy cascade may be observed even in a rough pipe, 
though it may be coexistant with an inverse cascade. Us- 
ing this method we can obtain the roughness dependence 



4 



of the friction factor in an enstrophy cascade dominated 
flow. The collapse of the friction factor curves using the 
enstrophy cascade variables is shown in Fig. The 
collapse is quite good, despite an apparent shallowness 
to the Blasius regime in the raw data. This shallow- 
ness is likely caused by the presence of a small amount 
of roughness, modifying the expected Re ' scaling at 
larger Reynolds numbers. We have neglected intermit- 
tency, which is negligible in 2D[2^. 

Relationship of the friction factor to the velocity profile:- 
FoUowing Prandtl 25|], we have calculated the mean ve- 
locity profile u{y) as a function of distance from a wall 
y, and for the enstrophy cascade this yields u{y) ~ j/" 
with a — 1/3, corresponding to the Blasius regime. For 
a general conserved quantity, a = (1 — (/))/(3 — 0). This 
relation depends on the zero roughness limit. In [2^ . 
it has been shown that roughness modifies the velocity 
profile so as to increase the apparent scaling exponent. 



Other work [271 123 ^^^o considers the influence of rough 
walls on the velocity profile and near-wall scaling. 

In our simulations of smooth-pipe enstrophy cascade 
turbulence, we have measured the velocity profile and 
found the power-law scaling exponent a = 0.323 ± 0.005 
between O.Oli? and O.li?, as shown in the inset of Fig. ([T]), 
close to the predicted a = 1/3. In the case of our rough- 
pipe simulations, the velocity profile yielded an exponent 
of 0.333 ± 0.002, significantly steeper than the predicted 
a = 1/7 that applies in the smooth, inverse cascade case. 
Our interpretation is that this is due to spectral contam- 
ination from an enstrophy cascade, as in the case of the 
simulations with random forcing that we presented. The 
momentum transfer theory integral has an upper limit 
that is comparable with the Kolmogorov lengthscale at 
low roughness, and so in that case the small-fc part of 
the energy spectrum controls the friction factor scaling. 
Because of this, we would expect to see a velocity profile 
consistent with the enstrophy cascade until the rough- 
ness or Reynolds number were high enough to place the 
crossover between the inverse cascade and contaminant 
enstrophy cascade below the scale of the roughness. 

Our results for the power-law Blasius regime in 2D 
enstrophy-dominated turbulence show convincingly that 
this regime is more than an empirical fit, and has a dy- 
namical significance. Our direct numerical simulations 
support the fundamental connection between spectral 
structure and friction factor scaling, which is manifested 
in the observed roughness-induced criticality. 
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supported by the National Science Foundation under 
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